Mon. Not. R. Astron. Soc. 000, [JHS] (2008) Printed 16 March 2011 (MN 1^1^ style file v2.2) 



The Tidal Downsizing hypothesis for planet formation and 
the composition of Solar System comets. 



o 



Sergei Nayakshin^'*, Seung-Hoon Cha^ and John Bridges^ 

^ Department of Physics & Astronomy, University of Leicester, Leicester, LEI 7RH, UK 

^ Space Research Centre, Department of Physics & Astronomy, University of Leicester, Leicester, LEI 7RH, UK 
* E-mail: Sergei.Nayakshin@astro.le.ac.uk 



Received 



6 



> 
in 
cn 
o 



ABSTRACT 

Comets are believed to be born in the outer Solar System where the temperature is 
assumed to have never exceeded T ^ 100 K. Surprisingly, observations and samples of 
cometary dust particles returned to Earth showed that they are in fact made of a mix 
of ices, as expected, but also of materials forged at high-temperatures (T ~ 1500 K). 
We propose a radically new view regarding the origin of the high-temperature pro- 
cessed materials in comets, based on the recent "Tidal Downsizing" (TD) hypothesis 
for planet formation. In the latter, the outer proto-planetary disc is gravitationally 
unstable and forms massive giant planet embryos (GEs). These hot (T ^ hundreds to 
2,000 K) and dense regions, immersed in the background cold and low density disc, 
are eventually disrupted. We propose that both planets and the high-T materials in 
comets are synthesised inside the GEs. Disruption of GEs separates planets and small 
solids as the latter are "frozen- in" into gas and are peeled off together with it. These 
small solids are then mixed with the ambient cold disc containing ices before being in- 
corporated into comets. Several predictions of this picture may be testable with future 
observations of the Solar System and exoplanets. 



1 INTRODUCTION 



Comets are icy bodies ^ km across that leave spectacular 
tails of material (dust) when ices are vaporised. Comets con- 
tain some of the most pristine materials from the dawn of the 
Solar System, and offer vital clues about its formation pro- 
cess. The composition of comets is confusingly diverse. Some 
of the materials found in cometary n uclei, e.g., amorphou s 
water and ammonia ices, have never (|Kawakita et al.ll2004l ) 
experienced temperatures above ~ 30 — 150 K, confirming 
their formation very far out, probably around the present 
day orbits of Uranus and Neptune. However, the mass frac- 
tion of crystalline silicates in the comae of the short period 
comet 81P/Wild 2, a nd in the ejecta of co met 9P/Tempel 1 
is t antalisingly high ( Zolenskv et al.|[2006l ). perhaps as high 
as (|Westphal et al.ll2009l ) ~ 0.5 - 0.65. This is surprising 
as some crystalline silicates such as o livine require temper - 
atures in excess of 1000 K to make (| Wooden et al.ll2007t ). 
although not all crystalline silicates form at high tempera- 
ture. 

In the "Core Accretion" (CA) paradigm of planet for - 
mation (|Safronovll 19721 : iPoUack et al ]|l996l : llda fc Linll2008l ) , 
the outer disc is a rather uninteresting and cold place where 
planet formation is not very like ly as the solid core forma- 
tion time scales are long ISafronov.. 196a : .Rafikov 20l3). The 
temperature outside _R 10 AU is generally expected to re- 
main below 100 K. Therefore, in the context of CA model, 
the presence of materials made at T ^ 1000 K in comets 
strongly sugg ests (|Wooden et al.ll2005l ) a radial transport of 



hight-T grains from the inner R ^ 1 AU regions into the 
outer , R ^ 10 — 30 AU regions. Detailed m odels of the pro- 
cess (|Gailll200ll : iHughes fc Armitage|[2O10l ) show numerous 
constraints necessary to satisfy in order to yield a significant 
enough outward transfer of solids. 

Here we sh ow that a set of recent ideas (|Bolev et al.l 
I2OIOI : I Navakshin 2010d Cibi). proposed as an alternative to 
the CA model for planet formation, as a by-product may 
naturally explain the otherwise puzzling composition of 
comets. The defining difference of the model from the CA 
scenario, as far as comet formation is concerned, is the non- 
unique radius-temperature relation in the disc. 

In the TD model, in a stark contrast to the CA pic- 
ture, the outer disc is the most important region for planet 
formation, as it is the birthplace of the giant planet em- 
bryos. These massive (~ 10 Jupiter masses) planet-forming 
gas clumps are very hot and dense not due to being close 
to the parent star or viscous disc heating, but simply due 
to contraction of the clumps. The clumps are in f act un- 
derm assive isolated "first cores" - embryos of stars (|Larsonl 
Il969h - that are not destined to develop into a lo w mass star 
due to the imposing presence of the parent star (|Navakshinl 
l2010bl ) that anchors the protostellar disc. It should thus not 
be surprising that these clumps manage to become as hot as 
~ 1000 K all on their own, at arbitrary distances from the 
parent star. 

Th e first cores (gas cl umps) are excellent sites for grain 
growth (|Navakshinll2010dl l3) and thermal processing of solid 
materials. Inside of these hot gaseous "ovens" , chemical com- 
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pounds can be baked into materials not normally expected 
to form at tens to hundreds of AU. Furthermore, a vital 
part of the TD model is the eventual disruption of gas em- 
bryos which release the planets back into the "ambient" disc. 
This disruption process, as we argue below, also releases the 
smaller thermally reprocessed solids back into the disc. The 
thermally reprocessed materials can then be rapidly mixed 
with the cold materials. As this is an in situ model, no out- 
ward transport of solids is required. 



2 THE TIDAL DOWNSIZING HYPOTHESIS 

The TD hypothesis is a new combination of earlier well 
known ideas and co ntains four impor tant stages (as illus- 
trated in Figure 1 of lNavakshinll2010al ): 

(1) Formation of gas clumps (which we also call giant planet 
embryos; GEs). A s the protoplanetary disc cannot fragment 
inside 7i ~ 50 AU l|Rafikovl2005l : lBolev et alll2006l '). GEs are 
formed at somewhat larger radii. The mass of t he clumps is 
estimated at Mqe ^ lOM j (10 Jupiter masses) (|Bolev et al.l 
I2OIOI : lNavakshinir2010dl ): they are intially fluffy and cool 
{T ~ 100 K), but contr act with time and become much 
hotter (|Navakshinll2010dl ). 



(2) Inward radial migration of the 
gravitational interactions with the 



clumps due to 
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(3) Grai n grow t h and sedimentation inside the c lumps 
llMcCreal Il960l: iMcCrea fc Williamsl Il965l : IBossI 1 19981 : 
IBoss et al.l 120021 1. If the clump temperature remains below 
1400 — 2000K, m assive terrestrial planet cores may form 
l|Navakshinll2010d '). with masses up to the total high Z ele- 
ment content of the clump (e.g., ~ 60 for a Solar met- 
alicity clump of total mass lOAfj). 

(4) A disru ption of GEs in the inner few AU due to 
tidal force s (iMcCrea fc Williamsl 1 19651 : iBolev et all l20ld : 
iNava kshiiJ '2010m or due to irradiation from the star 
l|Navakshin 2010b| ) can result in (a) a smallish solid core 
and a complete gas envelope removal - a terrestrial planet; 
(b) a massive solid core, with most of the gas removed - a 
Uranus-like planet; (c) a partial envelope removal leaves a 
gas giant planet like Jupiter or Saturn. For (b), an internal 
energy rel ease due to a massive core formation removes the 
envelope (jHandburv fc Williamslll975l : lNavakshinll2010j ). 

In contrast to the CA model, the TD scheme can- 
not work without a massive outer R ^ tens to a hun- 
dred AU region of the disc. The elements (3,4) from 
an earlier 1960s scenario fo r terrestrial planet formation 
llMcCrealll96d: iMcCrea fc WiUiams .196^ were rejected by 
iDonnison fc WilUamsl (| 19751 ) because step (1) is not possible 
in t he inner Solar System. Si milarly, the giant disc instabil- 
ity (|Kuipeij|l95ll : Bosslll998l') ca nnot operate at i? ~ 5 AU 
to make Jupiter ((Rafikovll2005l ). It is therefore the proper 
placement of step (1) into the outer reaches of the Solar 
System and then the introduction of the radial migration 
(ste p 2) that makes t his model physically viable. 

iNavakshinI ()201ll ) suggested that, as a bonus, the new 
hypothesis resolves an old mystery of the Solar System: the 



mainly coherent and prograde rotation of planets, which 
is unexpected in the CA framework since the planets are 
built by randomly oriente d impacts. Note, however, that 
iJohansen fc Lacerdal (|2010h shows that accretion of pebble- 
sized grains onto a planetary core could provide another 
explanation for the observed planetary spins. 

It is also not impossible (|Navakshinll2010bl ') that both 
the TD and the CA processes operate to sculpture the plan- 
etary systems we observe: the first in the early, gas-rich 
but s hort {t ^ 10^ yrs) embedded period (|Vorobvov fc Basul 
l2010h . and the second in the latter, much more quiescent 
phase t ^ a. few Myrs. In such a hybrid model the CA would 
kick-start with the benefit of the massive terrestrial cores 
pre-assembled in the early TD phase. 



3 RETAINING SMALL SOLIDS 

Although our arguments can be made c ompletely analytical, 
simulations of lCha fc NavakshinI (|2010l ) illustrate our model 
here. In the simulations, evolution of a massive 0.4 M0 gas 
disc around a 0.6 M© proto-star was followed for about 6000 
years. The massive disc becomes gravitationally unstable, 
develops spiral arms, which then fragment into clumps. The 
black solid curve in Figure [3^ shows the annuli-averaged 
and density weighted gas temperature from the simulation, 
defined as < pT > / < p > , whereas the solid curve in Figure 
^jp shows the corresponding density profile, < p >, as a 
function of radius R. The temperature and the density spikes 
correspond to the GEs in the simulation. To emphasise that 
even higher temperatures are present in the centres of the 
gas clumps, the red dashed curve in Figure [3^ shows the 
same as the black curve except for regions where density 
exceeds p > 10"^" g cm~"^, where the red curve shows the 
maximum temperature of the gas inside those regions. 

The black dashed curve in Figure ^jp shows the tidal 
density of the disc, pt = M,/{2TrR'^). Density of a disc 
marginally stable to the gravitational instability would fol- 
low the dashed curve. The "ambient" disc, i.e., the disc be- 
tween the gas clumps, has a density lower than pt and is also 
very cold, as expected. This confirms the two-phase division 
of the outer disc suggested above. 

A GE disruption should release gas and small solids 
with it back into the cold disc. Figure 2 shows the dust sed- 
imentation time scale as a function of dust particle radius, 
a, for GEs at three different ages from birth, 10^, lO-* and 
10^ y ears (dotted, solid and dash-triple-dotted, respectively) 
iNava kshin 2010bl). The grains are assumed to be located at 
ra, set to equal exactly the half radius of the embryo, Rge, 
from the embryo centre, but the results are almost indepen- 
dent of Vd- For example, the dashed curve is same as the 
solid one but calculated for = .02f?GE. A s the embryos 
are disrupted in ~ 10* to 10^ years (|Navakshin.2010bi' ) ■ small 
a <C 1 cm grains should be abundant in the gas envelope at 
the moment of disruption. Rapid radiative cooling and mix- 
ing with the cold background naturally deposits the high-T 
processed materials into the cold disc. 
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Figure 1. Solid curves: gas density (b) and temperature (a) av- 
er aged on annuli for numer ical simulation of a gas disc presented 
in lCha fc NavakshinI | |201C1) . The red curve showrs the maximum 
temperature found inside the clumps. The blue dotted curves 
show the corresponding temperature and density profiles in th e 
standard picture of planet formation llChiang fc GoldreichlllQQTl) . 
The nee d for radial transfer of solids by factors of ten or more is 
obvious l lGailll200ll : IWooden et al.ll2007t) in that theory. In con- 
trast, in the TD hypothesis, the outer disc has both hot and cold 
regions. Mixing of solids produced in these two components may 
yield a better explanation for the observed composition of comets. 



HOW ARE THE GAS CLUMPS DISRUPTED 
AT LARGE RADII? 



In the "bare-bone" version of the T D model, ste p 



disrupts t he env elope by tidal forces (iBolev et al.l 1201' 
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Navaks hinI l2010bl ) or due to stellar irradiation ( NavakshinI 
2010b|). Both of these effects are weak in the outer disc, e.s., 
beyond ^ 10 AU. What could disrupt the GE there, and is 
there any evidence for such disruption(s) in the Solar Sys- 
tem? 

As early as 35 years ago, iHandburv fc Williamsl (|l975l ) 
suggested that the massive core formation in Uranus and 
Neptune evaporated most of their hydrogen envelopes. The 
idea here is that the energy due to core formation is trapped 
inside the optically thick embryo, making it expand to sizes 
much larger than is expected in the simple analytical model 
of the gas clumps that do not take into accou nt the energy 
release by core formation (|Navakshinll2010bD . A more ex- 
tended gas envelope is then much easier to disrupt, even at 
tens of AU. 



To appreciate the argument, compare the binding en- 
ergy of the solid core with that of the GE. We expect the 
core of high-Z elements to have a density pc ~ a few g cm"''. 
The radial size of the solid core, Rcore ~ {3Mcore/4:TTpc)^^^ . 
The binding energy of the solid core is 
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Let us now look at the gas clump itself. The clump radius 
Rge ~ 0.8 AU at the age of t — 10* years, independently 
of its mass (iNavakshinlbOlObl ). Mge- Thus, the GE binding 
energy at that age is 
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The two are comparable for Mcore ~ lOMgj. Radiation hy- 
drodynamics simulations confirm such inter n al disr uption 
events: the run labelled MOaS in iNavakshi 3 (|2010d l made 
a 20 solid core that unbound all but 0.03 of the 
gaseous material of the original 10i\/j gas clump. 

If our model is right, then the outer ^ tens of AU Solar 
System must have produced at least one "naked" or almost 
so core as massive as 10 . There are actually two - Uranus 
and Neptune with core masses of ~ 13 and 15 M^, respec- 
tively. 

The self-disruption of GEs at tens or hundreds of AU is 
potentially observationally testable in exoplane tary systems. 
The Core Accretion model (|Pollack et al lliggel ) is unlikely to 
produce ^ 10 soli d cores that far out; the conve ntional 
disc instability model (|BossII1998| : iBolev et alll2006l ') would 
result in massive gas-dominated giant planets. Thus Super- 
Earth to Saturn mass planets on orbits of moderate eccen- 
tricity, if found at semi-major distances of ^ tens of AU, 
may be sign-posts of the gas clump self-destruction events. 

Another implication of our picture is that the composi- 
tion of comets and Neptune/Uranus may be related to some 
degree. The crystalline materials found in comets are mate- 
rials that did not contribute to the building of the gas giant 
planets. Estimates above show that forming solid cores is 
absolutely essential to the release of the high-T processed 
materials back into the cold disc. For the release to occur, 
the cores must be as massive as these outer planets, and so 
may have used up a significant fraction of the solids orig- 
inally present in the GEs. Therefore materials in comets 
that came from the same GE may be deficient in materi- 
als/elements abundant in Uranus and Neptune. 



5 IMPLICATIONS FOR THE ORIGIN OF 

OTHER HIGH-TEMPERATURE MINERALS 
IN THE SOLAR SYSTEM 

Other types of materials requiring high-temperature pro- 
cesses are chondrules and Calcium Aluminium-rich Inclu- 
sions (CAIs). Chondrules are igneous-textured, mm-size par- 
ticles, composed mainly of olivine and low-Ca py roxene set 
in a feldspathic or glassy matrix (e.g.. IScott fc Krot 200ij'). 
They are a major constituent of most chondrite groups (e.g., 
~ 80% of ordinary chondrites). The origin of chondrules is 
controversial, but in general they are believed to have formed 
as rapidly cooling molten silicate droplets. The maximum 
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temperatures are taken to be approximately around the liq- 
uidus temperatures of 1200-1500 K. The cooling rates re- 
main uncertain for the range of different textural types but 
it is thought that if chondrules had been molten for more 
than a few minutes they would not have preser ved the sort of 
volat ile abundances that they often contain (jYu fc HewinsI 
1 19981 ). 

CAIs are the light-coloured inclusions commonly found 
in carbonaceous chondrites. CAIs are more refractory- 
rich than chondrules. Their shapes are less regular, while 
common chondrules are more uniformly spherical. Radio- 
metric dating using the '^^Al — ^®Mg chronometer sug- 
;ests that chondrules star ted to form « 2 Ma after CAIs 
McKeegan fc Davis! 120051'). A Pb-Pb abso lute age for CAI 
formation is 4,567 Ma jAmelin et aLllioO^ ). 



Non-rotating embryo 



Thus CAIs are considered to predate chondrules. Nu- 
merous cho ndrule and C AI formatio n models include nebu - 
lar shocks (|Cassen|[l99^ . hghtn ing (iDesch fc Cuziill2000l ). 
jets fr om near the proto- Sun (|Shu et al.l l200ll ) and im- 
pacts (|Bridges et al.l Il998l ). No one model is universally 
accepted but the impact models ha ve the advantage o f 
producing abund a nt ch ondrules (e.g.. IScott fc Kro3l2005h . 
iHevev fc Sanderd (|2006l ) used the Ukely abundance of ^"^Al 
shortly after CAI formation in the early Solar System to 
show that nebular dust which rapidly accreted into ~ 60 km, 
or larger, planetesimals would start melting. Disruption of 
these planetesimals by impact would cause the sprays of 
melt droplets now seen preserved in chondrites. 

The GEs are pres ent only in the early "embedded" 
stage of star formation (|Vorobvov fc Basul l2010l : iNavakshinI 
l2010bl ). which is likely to last t ^ 10^ yrs. If this is true, and 
if the inferred age difference between the CAIs and chon- 
drules is real, then GEs are likely to be dispersed or become 
very massive giant planets by the time of chondrule forma- 
tion. 

If we assume that formation of CAIs is co-eval with 
the early GE-rich epoch of star and planet formation, then 
one may question whether GEs have anything to do with 
CAIs. We believe such a view is attractive because the tem- 
peratures near the solid core inside the gas embryos may 
l|Navakshinl l2010dl ) reach 1500-2000 K, e.g., high enough 
formation of CAIs. Vigorous convection ne ar the solid core 
jHelled fc Schubert! [ioog : lNavakshinll2010j ) probably drive 
strong shocks, which might be one way of produ cing CAIs. 
One-dimensional simulations o f iNavakshinI (|2010d ') also show 
melting/re- forming cycles for grains in some cases, e.g., see 
the right panel of his Fig. 8, the simulation M2a4. Physi- 
cally, the cycles result from a negative feedback loop. The 
accretion luminosity of the solid core increases as grains in- 
crease in size. However, this causes the inner GE regions to 
heat up, melting the grains. As grains become smaller, the 
rate of their accretion onto the solid core drops, and hence 
the luminosity of the core drops as well. The inner region 
cools down and the grains start growing again, repeating the 
cycle. Thus, in the TD model, this might explain the pres- 
ence of CAIs with original sizes up to ^ cm being found 
in comets. There is evidence for this from the Comet Wild2 
analyses and from IDPs (interplanetary dust particles). 
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Figure 2. Sedimentation time scales for grains versus their size, 
a. All tlie curves are calculated assuming tli e Mqe = lOMj gas 
embryo according to the model of Navaks hlnl l l2010bl '). but at three 
different embryo ages: 4 = 103, 10"' and 10^ years, for the dot- 
ted, the solid and the triple-dot-dashed curves, respectively. For 
these curves, the grains are located at half the gas clump ra- 
dius, whereas the dashed curve shows same as the solid curve 
but the grains are located at 0.02 the clump radius. The results 
are thus largely independent of the grain location inside the em- 
bryo. Grains smaller than a few mm to a few cm may remain 
suspended inside the embryo for a long time due to the gas drag 
forces. If the embryo is disrupted, these grains are released into 
the surrounding cold disc. 



6 CONCLUSIONS 

We have argued for an entirely different origin of the puz- 
zling comet compositions. Instead of assuming that hot 
T ~ 1500 K regions needed for thermal processing of hot 
minerals are located in the inner (_R ~ 1 AU) Solar Sys- 
tem, we identify them with the massive and appropriately 
hot gas clumps inside of which planets are born in the Tidal 
Downsizing hypothesis for planet formation. In the latter, 
the clumps are born at radii of many tens to hundreds of AU, 
and migrate inwards due to disc torques. The clumps are hot 
due to their self-gravity, and due to contraction caused by ra- 
diative cooling. We showed that small <g cm sized solids are 
suspended inside the gas clumps and ought to be released 
back into surrounding cold disc if the clump is disrupted. 
Disruption of the clumps in the outer Solar System requires 
a rapid formation of massive M > 10 M^ cores inside the 
gas clumps, which puffs up the gaseous envelope to the point 
of its removal. We tentatively identify Uranus and Neptune 
as two such cores that could have disrupted their gas em- 
bryos and contributed to building the comets in the Solar 
System. 

There may be chemical signatures of a casual link be- 
tween compositions of comets and the cores of the icy giants 
confirming (or rejecting) our model, perhaps testable with 
future results from the Rosetta mission. We also note that 
Tidal Downsizing hypothesis predicts massive solid cores 
(tens of Earth masses, e.g., planets like Uranus and Nep- 
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tune and possibly more massive) and Saturn-mass planets 
on semi-circular orbits at tens to hundreds of AU from the 
parent stars. Such planets are unlikely to be born in either 
the Core Accretion picture, where the core formation time 
is too I ons at _R ~ IQOAU, or in the disc instability model 
of fe.g.. lBossl[l99i l where the mass of th e fragment is muc h 
more likely to exceed that of Jupiter ('e.g.. lBolev et al]|201Cf ). 
This model-differentiating prediction may be testable with 
future observations of exoplanetary systems. 



ACKNOWLEDGMENTS 

The authors acknowledge illuminating discussions with and 
comments on the draft by Richard Alexander. Theoretical 
astrophysics research and cometary research at the Univer- 
sity of Leicester are supported by STFC rolling grants. 



REFERENCES 

Amelin Y., Krot A. N., Hutcheon L D., Ulyanov A. A., 

2002, Science, 297, 1678 
Boley A. C, Hayfield T., Mayer L., Durisen R. H., 2010, 

Icarus, 207, 509 
Boley A. C, Meji'a A. C, Durisen R. H., Cai K., Pickett 

M. K., D'Alessio R, 2006, ApJ, 651, 517 
Boss A. P., 1998, ApJ, 503, 923 

Boss A. P., Wetherill G. W., Haghighipour N., 2002, Icarus, 
156, 291 

Bridges J. C, Franchi I. A., Hutchison R., Sexton A. S., 
Pillinger C. T., 1998, Earth and Planetary Science Letters, 
155, 183 

Cassen P., 1996, in Chondrules and the Protoplanetary 
Disk, edited by R. Hewins, R. Jones, & E. Scott, 21-28 

Cha S.-H., Nayakshin S., 2010, submitted to MNRAS Let- 
ters 

Chiang E. I., Goldreich P., 1997, ApJ, 490, 368 

Desch S. J., Cuzzi J. N., 2000, ICARUS, 143, 87 

Donnison J. R., WiUiams I. P., 1975, MNRAS, 172, 257 

Gail H., 2001, A&A, 378, 192 

Goldreich P., Tremaine S., 1980, ApJ, 241, 425 

Handbury M. J., WiUiams I. P., 1975, AP&SS, 38, 29 

Helled R., Schubert G., 2008, Icarus, 198, 156 

Hevey P. J., Sanders I. S., 2006, Meteoritics and Planetary 

Science, 41, 95 
Hughes A. L. H., Armitage P. J., 2010, ApJ, 719, 1633 
Ida S., Lin D. N. C, 2008, ApJ, 685, 584 
Johansen A., Lacerda P., 2010, MNRAS, 404, 475 
Kawakita H., Watanabe J., Furusho R., et al., 2004, ApJ, 

601, 1152 

Kuiper G. P., 1951, in 50th Anniversary of the Yerkes Ob- 
servatory and Half a Century of Progress in Astrophysics, 
edited by J. A. Hynek, 357—1- 

Larson R. B., 1969, MNRAS, 145, 271 

Lin D. N. C, Bodenheimer P., Richardson D. C, 1996, 
Nature, 380, 606 

McCrea W. H., 1960, Royal Society of London Proceedings 
Series A, 256, 245 

McCrea W. H., Williams I. P., 1965, Royal Society of Lon- 
don Proceedings Series A, 287, 143 



McKeegan K. D., Davis A. M., 2005, Early Solar System 

Chronology, 431 — h, Elsevier B 
Nayakshin S., 2010a, ArXiv e-prints 
Nayakshin S., 2010b, MNRAS, 408, L36 
Nayakshin S., 2010c, ArXiv e-prints, 1007.4165 
Nayakshin S., 2010d, MNRAS, 408, 2381 
Nayakshin S., 2011, MNRAS, 410, LI 

Pollack J. B., Hubickyj O., Bodenheimer P., Lissauer J. J., 
Podolak M., Greenzweig Y., 1996, Icarus, 124, 62 

Rafikov R. R., 2005, ApJL, 621, L69 

Rafikov R. R., 2010, ArXiv e-prints 

Safronov V. S., 1969, Evoliutsiia doplanetnogo oblaka. 

Safronov V. S., 1972, Evolution of the protoplanetary cloud 
and formation of the earth and planets., Jerusalem (Is- 
rael): Israel Program for Scientific Translations, Keter 
Publishing House, 212 p. 

Scott E. R. D., Krot A. N., 2005, Chondrites and their 
Components, 143 — h, Elsevier B 

Shu F. H., Shang H., Gounelle M., Glassgold A. E., Lee T., 
2001, ApJ, 548, 1029 

Vorobyov E. I., Basu S., 2010, ArXiv e-prints 

Westphal A. J., Fakra S. C, Gainsforth Z., Marcus M. A., 
Ogliore R. C, Butterworth A. L., 2009, ApJ, 694, 18 

Wooden D., Desch S., Harker D., Gail H., Keller L., 2007, 
Protostars and Planets V, 815-833 

Wooden D. H., Harker D. E., Brearley A. J., 2005, in Chon- 
drites and the Protoplanetary Disk, edited by A. N. Krot, 
E. R. D. Scott, & B. Reipurth, vol. 341 of Astronomical 
Society of the Pacific Conference Series, 774 — h 

Yu Y., Hewins R. H., 1998, Geochim. Cosmochim. Acta, 
62, 159 

Zolensky M. E., Zega T. J., Yano H., et al., 2006, Science, 
314, 1735 



